A tensor instability in the Eddington inspired Born-Infeld Theory of Gravity 
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In this paper we consider an extension to Eddington's proposal for the gravitational action. We 
study tensor perturbations of a homogeneous and isotropic space-time in the Eddington regime, 
where modifications to Einstein gravity are strong. We find that the tensor mode is linearly unstable 
deep in the Eddington regime and discuss its cosmological implications. 
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I. INTRODUCTION 

An alternative theory of gravity was recently proposed 
in [H vi^hich attempted to extend Eddington's afHne the- 
ory of gravity to the matter sector ^] (see also 0, B])- 
The new theory is formulated in Palatini form in terms 
of the afRne connection P'^q,^ and a space-time metric, 
gap such that the gravitational action is given by 



\g^,u + KRf,y{T)\ - \y/g 



+S^[g,-^], 



where k = SttG, ^' denotes any additional matter fields. 
R^i, is the symmetric Ricci tensor constructed with P. 
This action can reproduce Eddington's original action at 
large values of kR and Einstein's at small values. 

It was shown in [l| that this theory has a novel be- 
haviour within sources. It was found that, in stars, com- 
pact objects and black holes, that the Eddington regime 
might not only lead to the avoidance of singularities but 
to significant modifications to the standard results of stel- 
lar astrophysics. The authors of ^ showed that it was 
possible to test the Eddington corrections to Newtonian 
gravity using Solar physics, while in Q it was shown 
that it should be possible to do the same around com- 
pact rotating sources. In Q it was shown that the mere 
existence of bound gravitational objects of a certain size 
led to stringent constraints on the free parameter of the 
theory while in Q it was shown that the theory could 
re-expressed as a bigravity theory, along the lines of Q . 

The Eddington regime will also arise in the very early 
Universe and it was shown in [l| that it may have led to 



a minimum scale factor. Hence, and more significantly 
than in black holes, the Eddington regime seems to pre- 
vent the formation of cosmological singularities, at least 
when seen in the context of homogeneous and isotropic 
space times. 

In this brief note, we look at one particular aspect 
of the cosmological behaviour of the Eddington regime 
by studying the structure and evolution of linear ten- 
sor mode perturbations, i.e. cosmological gravitational 
waves. We focus on tensor modes because they are purely 
gravitational, unlike the scalar modes that are a combi- 
nation of metric and density perturbations. This note is 
structured as follows: in section |lT] we lay out the main 
equations and derive the evolution equation for traceless, 
transverse perturbations (i.e. tensor modes) to a homo- 
geneous and isotropic space time; in section IIIII we find 
the solutions to the tensor equations in the Eddington 
regime and show that there is an instability; in section 
II VI we discuss the implication of our findings for this spe- 
cific model and the wider consequences of this tensor in- 
stability for bigravity theories and their equivalents. 



II. TENSOR MODES IN THE EDDINGTON 
BORN INFELD THEORY 

In this section we review some key features of \^ and 
introduce the procedure that will lead us to compute the 
tensor perturbations using the auxiliary metric q^v In 
Eq.([T]) we see that the metric and the connection are 
treated as independent variables, known as Palatini vari- 
ation. The resulting new set of evolution equations can 
be conveniently written in the form [l[ 
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{q-'r-Xg^'' - ~nT>^^{g,^), (1) 

qiiv - gt,u = nRt,u{q), (2) 



2 



where A = 1 + kA, and k is a constant with the inverse 
dimensions to that of the cosmological constant A. Here 
g^i^ is an auxihary rank-2 tensor, related to the original 
connection F via the Christoffel symbol. R^,^{q) is the 
associated curvature. The energy-momentum tensor is 
coupled to the metric These equations imply that 
T^'^ is conserved in the usual sense, T^^'^.^ — 0, where the 
covariant derivative is taken with respect to the metric 
g^i,. This property is not obvious from the equations of 
motion, but follows in a straightforward way from the La- 
grangian using invariance under general coordinate trans- 
formations. 

We can consider a perturbed homogeneous and 
isotropic space time by choosing the two metrics to be 
of the form 

g fj^udx^ dx'^ = —a^dff' + {Sij + hij) dx^dx^ , (3) 

qf^^dx^'dx" = -X^dr]'^ + -|- -f,j)dx'dx\ (4) 

where a, X and Y are solely function of conformal time, 
rj, both hij and ^ij are transverse and traceless, i.e. ha = 
-ii, = and diWi = 9,7'^ = 0. 

We can take the energy momentum tensor to be given 

by 

T^,^ ^ ip + P)Uf,u^ + Pg^^, (5) 

where u'^ = (1,0,0,0). The system can then be 

used to find the evolution equations at all times. We 
will be focusing on the very early Universe, when mat- 
ter could be described in terms of a relativistic perfect 
fluid. It was shown in fl!] that the modified Friedman- 
Robertson- Walker equation has the form: 



We obtain the field equations in the form: 



= - 



1 



Kp- 1 + 1)(3 - 

(1 



Kp){3 — kpY 



(3 + K2p2)2 



(6) 



where H = a'/a^ and for up « 1 we have ~ p/3. 
This equation has critical points for H{pb) = at a 
maximum density pB which depends on the sign of k. It 
is useful to write out the background equation that arise 
from the field equations: 



|A:y3| 


A 


P 




^— = 


\XY^\ 


A 


P 


Y^a-^ 





(7) 



To construct the perturbed field equations we need the 
following identities: 



\XY'' 

\XY'^ 
1 

P 



^0(7')), 
— Y^) where 



h'^) where h'^ = S^'S^'h 



„4 y27 +^2^ 



P 

K — h'^. 



a^Y^ 



But using Eq.(l7l) we have that 

XY^ 1 „■ A , ^ XY^ A 

which simplifies greatly to 

This is an intriguing result. Even though the tensor per- 
turbations in both the metric and the auxiliary metrics 
are multiplied by different conformal factors, they are 
identical in this theory. Furthermore, even in the Ein- 
stein regime, where X = Y = wc find that 7^ is non 
trivial and completely locked to the behaviour of hij . 

We can now proceed to construct the evolution equa- 
tion for hij to find 



Y' 

'y 

X^ 
Y2 



T 

-X^ 



hL 



f 1 



Y" X' Y' 

2 2 

Y X Y 



h,j = 0. 



(8) 



The spatial background field equation can now be used 
to find 



y2 = a2 



~Y ^ \Y) ^ YX 



which when replaced in Eq.® leads to our final expres- 
sion: 

2 



Y' X'\,, fx 

'y-xJ'-^^Hy 



Ph,j = 0. 



(9) 



This is a remarkably simple evolution equation for the 
tensor mode. In the Einstein limit it reduces to 



k^h, 



0, 



(10) 



as expected. 



III. EVOLUTION OF TENSOR MODES 

We now wish to see how this system evolves in the dif- 
ferent regimes. In the Einstein regime we find that the 
evolution is indistinguishable from Einstein gravity, even 
though the auxiliary metric is perturbed and present. 
The solutions to Eg. pU)) will be the standard ones, de- 
pending on the dominant source of energy momentum 
and can be found in any standard text book. For exam- 
ple in the radiation era, we find that 



oc 
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where 'HI'"'' (x) are Hankel functions of the mth kind. In 
the radiation era we have i/ = 1/2 and the solutions re- 
duce to a regular {sm{kri) / kr]) and irregular {cos{kr]) / krj) 
solution as krj — ^ 0. The familiar behaviour of gravita- 
tional waves fldl . emerges for very large kr] where 
we find decaying oscillatory solutions propagating at the 
speed of light. Crucial to this behaviour is that there is 
no time dependent factor multiplying the Laplacian (or 
Fourier space analogue, — fc^). 

It is in the Eddington regime that we find novel be- 
haviour and we will first focus on k > 0. In [l^ it was 
shown that the evolution of the background in the Ed- 
dington regime can be approximated by: 



a 

as 

2 




U = 



(11) 



where as is the minimum scale factor and t is physi- 
cal time. This primordial, non-singular behaviour was 
already alluded to in [1] and seems like an attractive, 
alternative explanation for early time. It resurrects an- 
other of Eddington's ideas of doing away with beginning 
and instead have the Universe indefinitely loitering in 
stasis in the distant past. We can see, however, from 
Eg. (1111) that the auxiliary metric does become singular 
as t — > — oo and it is this behaviour which will play a 
crucial role in the evolution of the tensor modes. 

To proceed we need to re-express our background 
quantities in terms of conformal time, rj which we can 
find by integrating dr/ = dt/a{t) using the above expres- 
sion to find: 

as 



1 — exp(Q!A77) ' 



where a = Ub •\/8/(3k) and Arj = rj — iji (where iji is a 
fiducial time) is negative (akin to when we write deSitter 
metric conformal coordinates we find that 77 < 0). We 
then have that 



Y2 



exp(Q!A77) 
1 — exp(Q;A77) 



= 2 



d,j ln(r V^) = drj In \fWJTj , 
exp(Q;A77) 



2a- 



1 — exp(Q!A77) ' 
The evolution equation for the tensor mode is 
exp(aA7y) ^, ^ exp(aA?7) 



2a 



[1 — exp(Q;A77)] 



[1 — exp(aA?7)] 



k^h, 



= 0, 
(12) 



A crucial difference with regards to tensor modes in Ein- 
stein gravity is that the pre factor of last term in Eq.(12) 
becomes singular. While in Einstein gravity that term is 
responsible for the acoustic, or wave-like behaviour of the 
tensor mode evolution, in Eddington gravity the Lapla- 
cian term is greatly suppressed. One possible interpreta- 
tion is that the speed of the gravitational wave goes to 
as i — —00. We can extract the asymptotic behaviour 
of the wave equation in the limit where A77 — > — 00 by 
discarding the last terms in Eq.(12) to find 



0, 



which we can solve to give hij oc Ai] + B. As announced, 
we have found an instability in the Eddington regime, in 
the asymptotic past. The consequences of such a result 
will be discussed in the conclusions of this report. 

We now turn to the case where k < where we found a 
bounce when the scale factor of the Universe reaches its 
minimum, non-singular value. It was shown in f]3\ that, 
if one choses a closed, positively curved spatial metric, it 
is possible to construct an oscillating (or Phoenix Uni- 
verse) which undergoes an indefinite number of cycles. 
Such a model should, in principle allow us to study the 
evolution of perturbations through the various cycles and 
shed light on some of the issues that have been raised in 
the study of cyclic cosmologies [H, [l3| • 

Again, we will work deep in the Eddington regime 
where the a, X and Y can be closely approximated by 
(assuming Euclidean geometry) 



3 k 




and we have assumed that the bounce occurs at i = 0. 
In conformal time we find that the scale factor can be 
expressed as 

a = asil + tan^(/3T?)], 



where /3 = aB\/2/ (3|k|). We then have 
^2^.4 1 



33/2 |tan(/3?7)|' 
^ |tan(/3?7)|. 



31/21 



(13) 



We can now take the Taylor expansion around 77 = to 
find that 



X2 
Y2 

Y' _ X' 



u 

V 



1 



2 
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The evolution equation for the tensor mode reduces to: 



2 P 



which can be solved with htj oc rj^ where 



(14) 



As in the case of k > we find an instability as a — ub, 
this time at the bounce; both solutions blow up as ?7~^/^ 
rendering such a space time unstable to tensor mode per- 
turbations. 



IV. CONCLUSIONS 

In this paper we found that, even though the back- 
ground evolution is resolutely non-singular, the overall 
evolution is still singular once one considers tensor per- 
turbations. It is an intriguing result and especially so in 
the case of k < where the tensor mode blows up at 
a finite time. The singular behaviour is clearly induced 
by the evolution of the homogeneous part of the auxil- 
iary metric (i.e. via X and Y) and is a completely novel 
effect, not present in conventional, Einstein gravity. 

A similar instability does arise in the radiation era, 
in Einstein gravity- as we saw above there is a decaying 
mode, proportional to l/krj at early times and clearly 
divergent as ?7 — 0. Yet, the current paradigm does 
not extend the radiation era all the way back to the Big 
Bang; there is an intervening period of de-Sitter expan- 
sion, known as cosmic Inflation, which itself has a finite 
duration. One can imagine invoking a more elaborate 
theory for Eddington theory in the asymptotic past but 
the simplest scenario, first proposed in [ij, does not work 
in its simplest incarnation. The instability is unavoidable 
for both K > and k < 0. 

How general is this behaviour? In and [1| , an alter- 
native view point for these theories of gravity was pro- 
posed, firmly placing them in the context of bigravity. 
Inevitably, in such theories, one finds more than one ten- 
sor mode which may or not be tightly coupled to each 



other, depending on whether the auxiliary metric has its 
own kinetic term. An analysis of perturbations in [TBj 
did not find such an instability in the scalar sector and 
this might be a hint that it is the particular form of the 
Eddington inspired Born Infeld theory that gives rise to 
such behaviour. Another possibility is the fact that we 
are considering a Palatini formulation of gravity where, 
as shown, pathologies occur that are absent in purely 
metric theories [lB| . We intend to look at general formu- 
lation of gravity theories [l3l to pin down the conditions 
in which such an instability arises. 

We may learn some lessons from the cosmological set- 
ting which may be applicable in other physical circum- 
stances. In particular, there has been some work on un- 
derstanding the process of gravitational collapse in Ed- 
dington gravity. The focus has been on spherical collapse, 
as one might expect. Clearly traceless, transverse modes 
may play an unexpected role and should be included if 
possible [l8| . Indeed by allowing for more general per- 
turbations it should also be possible trace the effect of 
nonlinear evolution of the tensor modes (coupled to ra- 
dial modes) to search if the singularity can be stabilised 
in the non-linear regime. 

Finally, there is of course the whole realm of gravita- 
tional waves to be explored in more general settings in 
these theories of modified gravity. Our analysis clearly 
hints at the possibility that interesting effects might arise 
in these theories in regions of density and curvature. This 
is clearly one of the new frontiers of modern gravitational 
physics which merits further exploration, specially given 
the developments in the study of inspiralling compact 
objects, pulsars and gravitational wave detection |19| . 



Acknowledgments 

We thank T.Baker and T. Clifton for useful discus- 
sions. This research is supported by STFC, Oxford Mar- 
tin School, BIPAC and the Basque Government through 
Ruth Lazkoz's research project AE-2010-1-31, Fundacion 
Pablo Garcia, FUNDEC, Mexico. MB was partially sup- 
ported by Fondecyt (Chile) Grants #1100282 and # 
1090753. 



[1] M. Banados, P. G. Ferreira, Phys. Rev. Lett. 105, 011101 

(2010). arXiv: 1006. 1769 [astro-ph.CO]] 
[2] A. Eddington, The Mathematical Theory of Relativity, 

CUP (1924). 

[3] S. Deser an d G. W. Gibbon s, Class. Quant. Grav. 15, 

L35 (1998) hep-th/9803049'. 
[4] D. N. Vollick, Phys. Rev. D 69, 064030 (2004) 
[5] J. Casanellas, P. Pani, I. Lopes, V. Cardoso, Astrophys. 

J., 747, 15, 2012 arXiv:1109.0249 astro-ph]] 
[6] P. Pani, V. Cardoso and T. Delsate, Phys. Rev. Lett. 107 



(2011) 031101 arXiv:1106.3569 [gr-qc]]. 

[7] P. Avelino, (2012) arXiv:1201.2544 

[8] T. Delsate, J. Steinhoff (2012) arXiv:1201.4989; 

[9] M. Banados, A. Gomberoff, D. C. Rodrigues and C. Sko- 

rdis, Phys. Rev. D 79, 063515 (2009) arXiv: 08 11.12701 

[gr-qc]]. 

[10] J. E. Lidsey, A. R. Liddle, E. W. Kolb, E. J. Copeland, 
T. Barreiro and M. Abney, Rev. Mod. Phys. 69, 373 
(1997) arXi v:astro-ph/9508078; . 

[11] A. R. Liddle, D. H. Lyth, Cambridge, UK: Univ. Pr. 



5 



(2000) 400 p. 

[12] J. Scargill, M. Banados, P. G. Ferreira (in preparation) 
2012. 

[13] R. H. Brandenberger, Phys. Rev. D80, 023535 (2009) 
[arXiv:09051514] 

[14] Y. S. Piao, Phys. Lett. B 677, 1 (2009) |arXiv:0901. 26441 
[gr-qc]]. 

[15] M. Banados, P. G. Ferreira and C. Skordis, Phys. Rev. 
D 79, 063511 (2009) ,arXiv:0811.1272 [astro-ph]]. 



[16] V. Vitaghano, T. Sotiriou, S. Liberati Phys. Rev. D82 

084007 (2010) arXiv: 1007.3937 
[17] T. Clifton, P. G. Ferreira, A. Padilla, C. Skordis Phys. 

Rep. 513 1 (2012) ,arXiv:1106.2 476 
[18] Y. -X. Liu, K. Yang, H. Guo and Y. Zhong, 

arXiv:1203.2349' 

[19] D. Psaltis (2008)' [arXiv: 0806. 1531] 



